Topological phase transitions with and without energy gap closing 
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Topological phase transitions in a three-dimensional (3D) topological insulator (TI) with an ex- 
change field of strength g are studied by calculating spin Chern numbers C ± (k z ) with momentum k z 
as a parameter. When \g\ exceeds a critical value g c , a transition of the 3D TI into a Weyl semimetal 
occurs, where two Weyl points appear as critical points separating k z regions with different first 
Chern numbers. For \g\ < g c , C ± (k z ) undergo a transition from ±1 to with increasing \k z \ to a 
critical value k z . Correspondingly, surface states exist for \k z \ < k z , and vanish for \k z \ > A;J. The 
transition at \k z \ = k z is acompanied by closing of spin spectrum gap rather than energy gap. 
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The quantum Hall (QH) effect Q, Q in a two- 
dimensional (2D) electron gas under a strong magnetic 
field provided the first example of topological state of 
matter in condensed matter physics, which cannot be 
described by the Landau theory of symmetry breaking. 
Thouless, Kohmoto, Nightingale, and Nijs (TKNN) re- 
vealed that the essential character of a QH insulator, dif- 
ferent from an ordinary insulator, is a topological invari- 
ant of occupied electron states [3[ or many-body wave- 
functions [J] . They related the Hall conductivity of the 
system to the first Chern number (or TKNN number), 
which is quantized when the Fermi level lies in an energy 
gap between Landau levels. In such systems, topological 
phase transitions can happen only by closing the energy 
gap. Gapless edge states must appear on the bound- 
ary between a QH insulator and an ordinary insulator, 
which is ensured by the topological invariant. Interest- 
ingly, Haldane Q proposed a spinless electron model on a 
2D honeycomb lattice with staggered magnetic fluxes to 
realize the topological QH effect without Landau levels. 

The quantum spin Hall (QSH) effect was first theoret- 
ically predicted by Kane and Mele [f| and by Bernevig 
and Zhang |7| , and then experimentally observed in HgTe 
quantum wells. [1,0] Unlike the QH systems, where time- 
reversal (TR) symmetry must be broken, the QSH sys- 
tems preserve the TR symmetry. The main ingredient is 
the existence of strong spin-orbit coupling, which acts as 
spin-dependent magnetic fluxes coupled to the electron 
momentum. The QSH state is characterized by a bulk 
band gap and gapless helical edge states on the sample 
boundary. [a-lllf The existence of the edge states is due 
to nontrivial topological properties of bulk energy bands. 
However, the bulk band topology of the QSH systems 
cannot be classified by the first Chern number, which al- 
ways vanishes. Instead, it is classified by new topological 
invariants, namely, the Zi index or the spin Chern 
numbers. |13l4l5| For TR-invariant systems, both Zi and 
spin Chern numbers were found to give an equivalent de- 
scription, [lj, ll5j The robustness of the Zi index relies 
on the presence of the TR symmetry. In contrast, the 



spin Chern numbers remain to be integer-quantized, in- 
dependent of any symmetry, as long as both the band 
gap and spin spectrum gap stay open. 1A\ They are also 
different from the first Chern number for the QH state, 
which is protected by the bulk energy gap alone. The 
spin Chern numbers have been employed to study the 
TR-symmetry-broken QSH effect. [16| 

The QSH system is an example of the 2D topological 
insulators (TIs). Its generalization to higher dimension 
led to the birth of 3D TIs. \v$2^ A 3D TI has a bulk 



band gap and surface states on the sample boundary. 
The metallic surface states provide a unique platform 
for realizing some exotic physical phenomena, such as 
Majorana fermions 22| and topological magnetoelectric 



effect. [23j, |24| The 3D TIs have been experimentally ob- 
served in Bii-^Sbz, Bi2Te3, and Bi2Se3 materials, [25- 
l3pj which greatly stimulates the research in this field. 
The 3D TIs with T R sy mmetry are usually classified by 
four Zi indices, [l?! [Hf and are divided into two general 
classes: strong and weak TIs, depending on the sum of 
the four Zi indices. In the presence of disorder, while 
the weak TIs are unstable, the strong TIs remain to be 
robust. The Zi indices are essentially defined only on 
the TR-symmetric planes in the Brillouin zone, and do 
not provide information about the distribution of surface 
states in the full momentum space. When the TR sym- 
metry is broken, the Zi indices become invalid. There- 
fore, a more general characterization scheme for the bulk 
band topology, which does not rely on any symmetry and 
can provide more information about the distribution of 
surface states, is highly desirable. 

In this Letter, for a 3D TI with an exchange field of 
strength <?, we consider a momentum component, e.g., 
k z , as a parameter, and analytically calculate spin Chern 
numbers C ± (k z ) for the effective 2D system. The phase 
diagram for C ± (k z ) obtained can describe the systematic 
evolution of the bulk band topology of the 3D TI with 
varying parameters, and provide more information about 
the surface states. For small g, C ± (fc z ) take values ±1 
for \k z \ < k z , and undergo a transition to at | k z | = k z . 
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Correspondingly, on a sample surface parallel to the z- 
axis, helical surface states exist in the region \k z \ < k°, 
and disappear for \k z \ > k°. At |fc z | = k°, the spin 
spectrum gap closes, but the energy gap remains open, 
which is distinct from a usual topological phase transi- 
tion, where the energy gap always collapses. When \g\ is 
greater than a critical value g c , a topological transition 
of the 3D TI into a Weyl semimetal occurs. Two Weyl 
points appear as critical points separating a QH phase 
of the effective 2D system for \k z \ < fc™ from an ordi- 
nary insulator phase for \k z \ > , indicating that their 
appearance is topological rather than accidental. Chiral 
surface states existing in the region \k z \ < fc™ give rise 
to the Fermi arcs. 

Let us start from the effective Hamiltonian proposed 
in Ref. 0: H = A 2 r x {k x a x + k y a y ) + M(k)r z + 
A\k z T x a z + <?<7 Z , which was used to describe the strong 
TI of Bi2Se3. Here, a m and r m (m = x,y or z) de- 
note the Pauli matrices in spin and orbital spaces, and 
M(k) = M - Bxk 2 z - B 2 k\ with k\ = k 2 x + k 2 y is the 
mass term expanded to the second order. In the last 
term, we include an exchange field of strength g, in order 
to study the TR-symmetry-broken effect on topological 
properties of the TI. For convenience, the momentum is 
set to be dimensionless, by properly choosing the units of 
parameters in the model, namely, M$, A\, A 2 , B\, and 
B 2 . 

Making a unitary transformation TL — W HU with 



U = M1 



i(l — r x )a z , we obtain 



H = A 2 (k x a x + k v a y ) + [M(k)r z + A x k z T x + g]a z . (1) 

The eigenstates of Hamiltonian (1) can be easily solved 
by first diagonalizing the operator in the square bracket. 
The four eigenenergies are obtained as 
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FIG. 1: Phase diagram determined from spin Chern numbers 
in the k z -g plane. The parameters are chosen to be Mo = 0.1, 
Ai = 0.7, A 2 = 1, Bi = 2.5, and B 2 = 14. 



remains to be a good quantum number. Correspondingly, 
nonzero C ± (k z ) indicate that edge states with the given 
k z must appear on the edge of the effective 2D system. 
The edge states at various k z essentially form surface 
states of the 3D sample. Therefore, the characteristics of 
the surface states can be determined from the calculation 
of the fc z -dependent spin Chern numbers C ± (fc 2 ). 

The spin Chern numbers for the effective 2D system 
are calculated in a standard way, which has been de- 
scribed in details in previous works. 15, By studying 
a special case of k z = 0, we find that the topological 
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properties of Eq. (1) can be described by the spin Chern 
numbers C {k z ) associated with r z . Here, r z corresponds 
to Ur z U' = r z a z in the original Hamiltonian H, and so 
can be considered as a spin operator, measuring the dif- 
ference of spin polarization between the two orbits. The 
eigenstates of projected spin operator Pt z P need to be 
calculated first, where P is the projection operator into 
the valence bands. Since Pt z P commutes with momen- 
tum operator, its eigenstates can be obtained at each 
momentum k separately. The eigenvalues of Pt z P are 
given by 



where A(k) = \J M 2 (k) + A 2 k 2 , and subscripts ± indi- 
cate two valence (conduction) bands with superscript v 
(c). The electron wavefunctions in the valence bands are 
given by 



P±(k) = 0±(k)®x±(k) . 



(3) 



Here, </>+(k) = [sgn(Aifc z ) cosafe, sinafc] T and </>-(k) = 
[sgn(j4iA; z ) sina^, — cosak] T are wavefunctions in the t z 
space, and x± — ( kx k ^ ky sin 9^ , — cos 0^) T are wavefunc- 
tions in the a z space, with 2a k = ctg" 1 [M(k)/|Aifc z |] 
and 20± = ctg-^Cg ± A(k))/|A 2 fcx|]. 

The basic idea of our theoretical calculation is ex- 
plained as follows. We consider one of the momentum 
components, e.g., k z as a parameter. For a given k z , Eq. 
(1) is equivalent to a 2D system, for which spin Chern 
numbers C^(k z ) can be defined. For a semi- infinite sam- 
ple of the 3D TI with its surface parallel to the z axis, k z 



f*(k) = ±-\/cos 2 2a k + cos 2 (#+ - 9-) sin 2 2a k . (4) 

The corresponding eigenfunctions are denoted by 4' ± (k), 
whose expressions are lengthy and will not be written 
out here. The spin Chern numbers are just the Chern 
numbers of the two spin sectors formed by 4' ± (k), i.e., 
C^fc,) = £jdk x dk y e z ■ [V a x (* ± (k)|V 2 |* ± (k)], 
where V 2 is the 2D Laplace operator acting on (k x , k y ). 
By some algebra, C ± (fc z ) are derived to be 



C ± {k z ) = ±- 



5gn(B 2 ) + sgn(g(fc,)sgn[P(fc,)] ± 5 )] (5) 



with P(k z ) = M - B x k\ and Q(k z ) = ^P 2 (k z ) + A\k 2 z . 

Equations (2), (4), and (5) are the main analytical 
results of this work. For g = 0, the spin Chern numbers 
at k z = reduces to 6^(0) = ±|[sgn(-B 2 ) + sgn(M )]. 
(0) are nonzero when B 2 and Mq have the same sign, 
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FIG. 2: Energy gap Ae (a) and spin spectrum gap A r (b) 
as functions of momentum k z for different values of g. The 
other parameters are the same as in Fig. 1. 



and vanish otherwise. C ± (0) play a role similar to the Z 2 
index. Nonzero C (0) ensure that surface states exist in 
the vicinity of k z = on a surface parallel to the z axis. 
Without loss of generality, we will focus on the parameter 
region of B 2 > and M > 0, to which Bi 2 Se3 belongs. 
We wish to emphasize here that when k z is considered as 
a parameter, the effective 2D Hamiltonian (1) breaks the 
TR symmetry for any k z ^ 0, even if y = 0, as its TR 
couterpart is at — k z . Therefore, while the ^-dependent 
spin Chern numbers given by Eq. (5) remain to be valid 
at any k Zl a Z 2 index cannot be defined for any k z 0. 
Eq. (5) allows us to extract more information about the 
basic characteristics of the surface states. 

A typical phase diagram for the spin Chern numbers 
in the k z versus g plane, as determined by Eq. (5), is 
plotted in Fig. 1. For simplicity, A 2 is taken to be the 
unit of energy. For small \g\, C (k z ) — ±1 at small 
\k z \, corresponding to a QSH phase of the effective 2D 
system, and drop to with increasing \k z \ to a critical 
value k z = M jB\ , as indicated by the dotted lines. 
The system becomes an ordinary insulator for \k z \ > k z . 
When \g\ is greater than a critical value g c , the effective 
2D system enters a QH phase with a nonzero total (first) 
Chern number C(k z ) = C+{k z ) + C~{k z ) = 1 if g > 0, 
and —I if g < 0. The boundary enclosing this phase 
is determined by equation g = ±Q(k z ), as indicated by 
the solid curves. The critical exchange field is given by 
g c = min[Q(fc z )], and g c ~ 0.1 for the parameters used 
in Fig. 1. 

It is interesting to see how the energy gap A^fcj) = 
min{E±(k) — £/±(k)]|fc a and the spin spectrum gap 
A T (k z ) = min[^ + (k) — £~(k)]|fc z behave on the boundary 
(solid and dashed lines) between different phases. From 
Eqs. (2) and (4), we find that on the dotted lines in Fig. 
1, the spin spectrum gap closes at k x = k y = 0, but the 
energy gap remains open. On the contrary, on the solid 
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FIG. 3: Profiles of energy spectrum for a semi-finite sample 
of the 3D TI for four different k z at g = 0.05. The other 
parameters are the same as in Fig. 1. 
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FIG. 4: (a) Maximum level spacing SE between the edge 
states and bulk states as a function of k z /k z . (b) A schematic 
of the distribution of surface states in 2D momentum space. 



boundary lines, the energy gap closes at k x = k y = 0, 
but the spin spectrum gap remains open. and A T as 
functions of k z for several values of g are plotted in Figs. 
2(a) and 2(b), respectively. We notice that at g = 0.14 
and k z = k z = 0.2, the energy gap and spin spectrum 
gap vanish simultaneously. This is because the dotted 
and solid boundary lines in Fig. 1 intersect just at that 
point. 

To study the surface states directly, we construct a 
tight binding model on a cubic lattice with two spins and 
two orbits on each site, which recovers the Hamiltonian 
Eq. (1) in the continuum limit. A semi-infinite sample 
with its surface parallel to the y-z plane is considered, 
where k y and k z remain to be good quantum numbers. 
The calculated energy spectrum for g = 0.05 is plotted 
as a function of k y for four different values of k z in Fig. 
3. Although g ^ breaks the TR symmetry, the surface 
states remain to be gapless at k z = 0, because t z in Eq. 
(1) is conserved at k z = 0. From Fig. 3, it is found that 
for k z = 0, k z /3, and 2k z /3, surface states always exist 
in the bulk energy gap, but no surface states appear at 
k z = k z . To see the evolution of surface states with 
k z more clearly, we define a maximum level spacing SE 
between the surface states and bulk states, as illustrated 
in Fig. 3. In Fig. 4(a), SE is plotted as a function of k z . 
One can see that SE decreases with increasing k z , and 
drops to nearly at k z > k z . Therefore, we conclude 
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FIG. 5: Profiles of energy spectrum of a semi-finite sample 
for four different k z at g = 0.2. The other parameters are the 
same as in Fig. 1. 



that surface states exist only in the region \k z \ < k°, an d 
vanish for \k z \ > k°, which is well consistent with the 
phase diagram of Fig. 1. 

The critical momentum marks a topological phase 
transition, characterized by the change in the spin Chern 
numbers, with disappearance of surface states as an ob- 
servable consequence. This topological phase transition 
is unconventional, which is accompanied by closing the 
spin spectrum gap rather than the energy gap. While 
the closing of the spin spectrum gap may not be ob- 
served directly, we find that the average of operator r z 
in H, namely, 52p=±(<PpQ s )\ UW z U\(fp(k)) oc cos(2a fc ), 
changes its sign at k = (0, 0, ±fc°), which indicates a re- 
versal of the orbital polarization, leading to a bulk physi- 
cal observable at the transition. For the strong TI under 
consideration, the nontrivial bulk band topology can be 
examined in any direction. For example, by considering 
k y as a parameter, we can calculate the spin Chern num- 



bers in the effective 2D space of (k x ,k z ), and obtain a 
critical ky = y/Mo/Eh. Combining ky and k z together, 
we can depict the region in the k y -k z plane, where topo- 
logical surface states exist on a surface parallel to the y-z 
plane, as shown in Fig. 4(b). 

For \g\ > g c , the 3D system enters another topologi- 
cal phase characterized by a nonzero total Chern num- 
ber C(k z ) for small JfcJ, which is essentially a Weyl 
semimetal phase. [31 - 34 1 The quantum phase transition 



from the 3D TI to the Weyl semimetal with tunning g 
can be understood as a topological transition, at which 
the energy gap closing causes one of spin Chern numbers 
C ± (fc z ) to vanish, while the other remains to be quan- 
tized. At a given g, from equation g = ±Q(k z ) for the 
solid boundary lines in Fig. 1, one can obtain two crit- 
ical values of k z , namely, ±/c", as indicated in Fig. 1. 
k = (0, 0, ±k™) are a pair of Weyl points, or 3D Dirac 
points, at which the conduction and valence bands touch. 
The appearance of Weyl points is usually attributed to 
accidental degeneracy or symmetry. 3lM34| The obtained 



system for | k z \ < k™ from an ordinary insulator state for 
\k z \ > at which the energy gap must vanish. Chang- 
ing parameters cannot open energy gaps at the two Weyl 
points, unless they come together, so that the QH state in 
between is destroyed. At g — 0.2, the calculated energy 
spectra for four different k z as functions of k y are plotted 
in Fig. 5 for a semi-infinite sample with its surface paral- 
lel to the y-z plane. It is found that for k z < fc™, chiral 
surface states appear in the energy gap; the conduction 
and valence bands touch at k z = k™; and for k z > A;™, 
the energy gap reopens, but surface states no longer ex- 
ist. These results are in good agreement with the phase 
diagram in Fig. 1. The chiral surface states appearing in 



the region of \k z \ < k™ give rise to the Fermi arcs. [31 
l34l | The Weyl semimetals with Fermi arcs have been paid 
much attention recently. Our work suggests that a pos- 
sible route to realize such an interesting 3D topological 
state of matter is through magnetic doping in 3D TIs. 
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phase diagram in Fig. 1 sheds more light on its topolog- 
ical origin. At a given £/, the Weyl points appear as two 
critical points separating a QH state of the effective 2D 
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